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Nb nano superconducting quantum interference devices with high spin sensitivity for 

operation in magnetic fields up to 0.5 T 
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We investigate electric transport and noise properties of microstrip-type submicron direct current 
superconducting quantum interference devices (dc SQUIDs) based on Nb thin films and overdamped 
Josephson junctions with a HfTi barrier. The SQUIDs were designed for optimal spin sensitivity 
Sp upon operation in intermediate magnetic fields B (tens of mT), applied perpendicular to the 
substrate plane. Our so far best SQUID can be continuously operated in fields up to B m ±50 mT 
with rms flux noise Sj w < 250n<I > o/Hz 1 ' 2 in the white noise regime and spin sensitivity Sj < 



29/xb/Hz 1 ' 2 . Furthermore, we demonstrate operation in B — 0.5 T with high sensitivity in flux 

~ 79/^b/Hz 1 ' 2 . We discuss strategies to further 



Sl/l w 680n$ /Hz 1/2 and in electron spin S, 



improve the nanoSQUID performance. 



PACS numbers: 85.25.CP, 85.25.Dq, 74.78. Na, 74.25.F- 74.40.Dc 



Recent developments in miniaturized submicron-sizcd 
direct current (dc) superconducting quantum interfer- 
ence devices (SQUIDs) are motivated by the need of sen- 
sitive detectors for small spin systems such as molecu- 
lar magnets^— and magnetic nanoparticlcs,— cold atom 
clouds^ or single electrons and atoms^ and improved res- 
olution in scanning SQUID microscopy^— As a common 
approach, nanoSQUIDs based on constriction Joseph- 
son junctions (JJs) have been used^^r— achieving root 
mean square (rms) flux noise power S$ down to a 
few lOOn^o/Hz 1 / 2 ($0 is the magnetic flux quantum) 
in magnetically shielded environment.— However, con- 
striction JJs, even if resistively shunted, often show hys- 
teretic current- voltage-characteristics (IVCs). This ham- 
pers continuous SQUID operation as required for the in- 
vestigation of magnetization dynamics of magnetic par- 
ticles and the use of common SQUID electronics, de- 
veloped for readout of very sensitive dc SQUIDs with 
nonhysteretic JJs. Furthermore, the noise properties of 
constriction JJs are not well understood, which luakes 
SQUID optimization difficult. 

An alternative approach is the use of submicron 
superconductor-normal conductor-superconductor (SNS) 
sandwich-type JJs, which offer large critical current den- 
sities in the 10 5 A/cm 2 range and which are intrinsically 
shunted, providing nonhysteretic IVCs without the need 
of bulky external shunt resistors J^ In a standard thin film 
SQUID geometry, the SQUID loop and the JJ barrier are 
in the plane of the thin films. For detection of magneti- 
zation reversal of a small magnetic particle, one applies 
an external magnetic field in the plane of the SQUID 
loop and detects the change of the stray field coupled to 
the SQUID upon magnetization reversal, without cou- 
pling the external field to the SQUID. However, in this 
case, the applied field also couples magnetic flux into the 



JJ barrier and reduces its critical current, which in turn 
reduces the SQUID sensitivity. In order to avoid this 
problem, in this letter we present results on a modified 
SQUID design, which takes advantage of the multilayer 
technology used for SNS JJ fabrication. This approach 
allows for a further reduction of the SQUID inductance 
and hence improved SQUID sensitivity and at the same 
time operation in higher magnetic fields. 

The Nb thin film dc SQUIDs have a microstrip geom- 
etry, i.e. the two 250 nm wide arms of the SQUID loop 
lie directly on top of each other. The 200 nm thick bot- 
tom and 160 nm thick top Nb layers are separated by 
a 225 nm thick insulating Si02 layer and are connected 
via two JJs with areas 200 x 200 nm 2 and a nominally 
24 nm thick HfTi barrier (see Fig. [I]). HfTi was chosen 
as a barrier material as, among other binary materials, 
it provides a relatively high resistivity, does not become 
superconducting at 4.2 K and is compatible with our fab- 
rication technology. For details on sample fabrication 
and J J properties we refer to [19I - 421J J . The size of the 
SQUID loop is defined by the 1.6 //m spacing between 
the JJs and by the SiC>2 interlayer thickness. In contrast 
to earlier work^ for this geometry a sufficiently large 
magnetic field B can be applied perpendicular to the sub- 
strate plane without inducing a significant magnetic flux 
penetrating either the SQUID loop or the junction bar- 
rier. Furthermore, this design provides a very small area 
of the SQUID loop and hence a very small SQUID induc- 
tance L of a few pH or even lower. This is essential for 
reaching ultra-low values for the spectral density of flux 
noise power S&2& 

For current and flux biasing, additional 250 nm wide 
Nb lines are connecting the SQUID in a cross-shape ge- 
ometry, and a bias current lb, flowing from the top Nb 
layer through the JJs to the bottom Nb layer, can be 




FIG. 1: (Color online) Scanning electron microscopy (SEM) 
image of SQUID 2. Open (yellow) squares indicate positions 
of JJs. Arrows indicate current paths for bias current 7b 
(dashed: symmetric bias; solid: asymmetric bias) and modu- 
lation current 7 mo d (dotted). 



applied either in a symmetric or asymmetric configura- 
tion; see Fig. [TJ For simplified readout we use asym- 
metric current bias in the following. A magnetic flux $ 
can be coupled into the SQUID loop by applying a mod- 
ulation current 7 mo d across the bottom Nb layer ("flux 
bias line"). This enables flux biasing the SQUIDs at the 
optimum working point without the need of an external 
coil. Furthermore, the flux bias line can also be used to 
provide a feedback flux for SQUID operation in a flux 
locked loop. However, in this work, the SQUIDs were 
always read out open loop. 

We investigated various SQUIDs which were fabri- 
cated in two different runs on separate wafers. Below 
we present results for two devices, SQUID 1 from wafer 1 
and SQUID 2 from wafer 2. The main difference in the 
design of these devices is the different lengths ~ 2.5 /xm 
(SQUID 1) and - 5/jm (SQUID 2) of the narrow bias 
lines, running from the center of the SQUID to the 
4 /jm wide connection lines further away from the SQUID 
(cf. Fig. []}. This variation has a strong impact on the 
SQUID performance in applied magnetic fields, as will 
be shown below. 

All data were taken at temperature T = 4.2 K. We 
first present results of transport and noise measurements 
of the two SQUIDs in an electrically and magnetically 
shielded environment. Since both devices showed quali- 
tatively the same behavior, we only give a detailed anal- 
ysis of SQUID 2 and summarize the main parameters ex- 
tracted for both devices in Table |U Regarding absolute 
values, a major difference between both devices are the 
values for maximum critical current I c and normal re- 
sistance i?N, which probably is due to variations in the 
HfTi barrier thickness for the devices fabricated in differ- 
ent runs. Devices fabricated from the same run showed 
a spread in I c and i?N values of ±10 %. 

Figure l^a) shows IVCs of SQUID 2 for $/$ = 0, 1/4 
and 1/2. The IVCs are nonhysteretic with I c = 227 /i A 



and i?N = 250 nrfi, yielding a characteristic voltage 
V c = 7 c i? N = 57 (j,V. The IVC at $/$ = 1/2 exhibits a 
small bump for low voltages. This bump appears in all 
our devices and is presumably a property of the quasipar- 
ticle current rather than a LC resonance of the SQUID. 
The inset of Fig. [5Ja) shows the modulation I c (Imod) 
for positive and negative bias current. From the mod- 
ulation period we obtain the inverse mutual inductance 
AI-^ 1 = 2.73mA/$o- From the modulation depth we 
find a screening parameter /3l = 2IqL/$>o = 0.25. By 
assuming that both JJs are identical, i.e. I c = 2Iq, we 
determine the SQUID inductance L = 2.3 pH. 

The T^(imod) modulation for different bias currents, 
plotted in the inset of Fig. [2Kb), yields a maximum trans- 
fer function V 9 = dV/d$ = 164 ^V/$ for h = 230 /xA. 
The shift in I c {I mo d) and V(-fmod) for positive and neg- 
ative bias currents can be attributed to the asymmetric 
current bias, which leads to an inductance asymmetry 
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FIG. 2: (Color online) Transport and noise characteristics of 
SQUID 2. (a) IVCs for different flux <!>; inset shows measure- 
ment (solid black lines) and numerical simulation (dotted red 
lines) of 7c (/mod), (b) Solid black line: Spectral density of rms 

1 /2 

flux noise S$ (/) at optimum working point (lb — 230 ^iA, 
/mod = 243 /J, A). Dotted (red) line: fitted spectrum; dashed 
(red) line indicates white noise level for fitted curve. Inset 
shows V(J mod ) for I b = ±(150 . . . 300) /xA (in lOixA steps). 



TABLE I: Parameters of SQUID 1 and SQUID 2. 



I c i?N ^c^N PL 

( M A) (mfl) (/iV) 



L M 



-i 



V* 5, 



1/2 



<E>,w 

( P h) (g^) (g) (A; 



SQUID 1 129 385 50 0.19 3.0 2.63 154 260 
SQUID 2 227 250 57 0.25 2.3 2.73 164 200 



ckl = {L2 — Li)/(Li + L 2 ); here L 1 and L 2 are the induc- 
tances of the two SQUID arms. The measured I c (^mod)- 
characteristics are fitted well by numerical simulations 
based on coupled Langevin equations^ with a noise pa- 
rameter T = 2irk B T/I $ = 1.55 ■ 10" 3 (fc B is the Boltz- 
mann constant) and «l = —0.35 (see inset of Fig. [2a) ; 
dotted lines). 

Using a commercial SQUID amplifier with a voltage 
noise S-J ~ iOpV/Hz 1 / 2 and a -3dB cutoff frequency 



/ c ~ 30 kHz, we measured the spectral density of the rms 
flux noise S$ (/) = S v (/)/ |V$| at the optimum work- 
ing point (see solid line in Fig. Hfb))- Here the SQUID 
amplifier contribution was subtracted. We observe a sig- 
nificant low-frequency excess noise, which wie assign to 
Iq fluctuations in the JJs. Since the low-frequency ex- 
cess noise extends to well above 1 kHz and due to the 
limited bandwidth of the SQUID amplifier, we do not 
see a clear white noise region in the spectrum. By fit- 
ting the experimental data (dotted line in Fig. [2](b)), we 
derive a low-frequency noise contribution S^ f oc l// a 
with a = 0.5 and sffif = 1Hz) = 3.7 ^o/Hz 1 / 2 and a 

white noise contribution Sj w = 200 n$ /Hz 1,/2 (dashed 
line inFig.IHb)). 

1/2 

In order to determine the spin sensitivity Sp = 

5 $ ' /4>fj, of our SQUIDs, we calculated the coupling fac- 
tor (j)^, using a routine based on the numerical solution of 
the London equations for the given SQUID geometry^ 
Here, (f>^ = $//x is the magnetic flux $ per magnetic 
moment \jl\ = [i coupled by a magnetic particle to the 
SQUID loop. Very recently, the validity of this approach 
has been verified experimentally by measuring the mag- 
netic coupling of a Ni nanotube to a Nb nanoSQUID 
which had the same geometry as SQUID 2 2^ For a point- 
like magnetic particle with fl perpendicular to the sub- 
strate plane, placed at a lateral distance of 10 nm from 
the lower edge of the upper Nb SQUID arm at the center 
of the loop, we obtain 4>^ = 8.6n$o//XB (/■'b is the Bohr 
magneton). Along with the obtained value of the rms 



flux noise S, 



1/2 



1/2 

sensitivity to S^ 



200 n^o/Hz 1 / 2 we calculate the spin 

= 23^b/Hz 1 / 2 . 
To investigate the SQUID performance in a magnetic 
field B applied perpendicular to the substrate plane we 
mounted SQUID 1 on a high-precision alignment system 
(one rotator, two goniometers). B is generated by a 
superconducting split coil running in persistent mode 



to suppress field noised Figure [3ja) shows I C (B) for 
SQUID 1 after the alignment process for a field sweep 
sequence as indicated by labels 0-6. The observed hys- 
teresis for \B\ < 45 mT is ascribed to entry and trapping 
of Abrikosov vortices in the 4/xm wide connection lines, 
cf. inset of Fig. [3Jb). The steep jump in I c at B fa 45 mT 
can be assigned to a vortex entering the narrow Nb leads 
very close to the SQUID loop, as confirmed recently by 
magnetic force microscopy on a similar Nb nanoSQUID 
(with layout of SQUID 2)2^ Subsequently, we reduced 
the linewidth of the connection lines of SQUID 1 from 
4/xm to ~ 500 nm by focused ion beam (FIB) milling^ 
see inset of Fig. [3jb). For the repatterned device, the 
maximum I c was reduced by ~ 10%, probably due to a 
slight degradation of the JJs during FIB milling. More 
importantly, I c became almost independent of B, and 
within B w ±50 mT the magnetic hysteresis disappeared, 
cf. Fig. [3Jb). At B w 50 mT we still observed the jump 
in I c due to vortex entry in the narrow Nb line close 
to the SQUID. This indicates that the linewidth of the 
Nb wiring close to the SQUID may limit the range of 
operation to \B\ < 50 mT. However, as will be shown 
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FIG. 3: (Color online) I C (B) data of SQUID 1 for field sweep 
sequence 0-6 (a) and 1-3 (b) after removing Nb areas by FIB 
milling as indicated by hatched (yellow) rectangles in the inset 
(SEM image). 



below, even after vortex entry, by proper realignment of 
the applied magnetic field direction, which compensates 
the stray magnetic flux induced by trapped vortices, I c 
can be restored and low flux noise can be retained. 

We now turn to SQUID 2, which has much longer nar- 
row bias lines. Figure SJa) shows I C (B) for a field sweep 
46 mT ->■ -46 mT ->■ 55 mT (1-3). Again I c is almost 
independent of B for \B\ < 50 mT and, as before, we 
find a jump in I c at B rs 50 mT due to a vortex enter- 
ing the narrow bias lines. The vortex can be removed by 
sweeping back the field as indicated by the curve (3-4) 
in Fig. gja). 

For SQUID 2 we performed noise measurements as de- 

1/2 

scribed above to determine S§ w at several values of B 
from to 50 mT, without any jump in I c , see inset of 
Fig. Ha). For B = 0, £&£ w 220n$ /Hz 1 / 2 , which 
is slightly higher than the value obtained in the low- 
field setup. We attribute this to external disturbances 
from the unshielded environment in the high-field setup 



(cf. noise spectrum in Fig.[4|b), black line). As indicated 
in the inset of Fig. |4|a), the white noise level increases 
only slightly with B to Sl(l « 250n$ /Hz 1 / 2 at B = 
50 mT (cf. noise spectrum in Fig.HJb)), still correspond- 
ing to a very small spin sensitivity Sp/ ps 29//b/Hz 1//2 
(in the white noise regime) . We assign this behavior to a 
minor decrease of I c due to an imperfect alignment of the 
device relative to B. At B — 55 mT, i.e. after the jump in 
I c occurred and after realigning the SQUID by maximiz- 
ing J c , we obtain a similar value Sj « 240 n^o/Hz 1 / 2 as 
for B = 50 mT. Following the same procedure of realign- 
ment, we were able to operate the SQUID in magnetic 
fields up to B = 0.5 T, yielding the noise spectrum as 



,1/2 



,1/2 



shown in Fig. BJb), with S^^ w 680n$o/Hz 1/ ', corre- 
sponding to SfJ « 79/zb/Hz 1 / 2 . Note that all spectra 
feature excess low-frequency noise peaks, which are pre- 
sumably due to mechanical vibrations of the setup. 
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In conclusion, we fabricated and investigated Nb 
nanoSQUIDs based on a trilayer geometry which were 
optimized for stable operation in comparatively large 
magnetic fields. Very low white flux noise values down 
to Sj w ~ 200 nf&o/Hz 1 / 2 have been achieved in a 

1/2 

shielded environment yielding a spin sensitivity Sy w 
23/j.b/Hz 1 / 2 . Concerning the suitability to applied mag- 
netic fields, we successfully redesigned the layout of 
SQUID 1 via FIB milling and implemented these findings 
into the design of SQUID 2. We demonstrated stable op- 
eration in a field range of B « ±50 mT with a marginal 
increase in white flux noise and spin sensitivity with B 
(Sl{l < 250n$ /Hz 1 / 2 and S l J 2 < 29 /is/Hz 1 / 2 ). More- 
over it was shown that SQUID 2 can maintain high sen- 

1/2 

sitivity in large fields up to B = 0.5 T with S$ « 
680 mJo/Hz 1 / 2 and S^ 2 w 79^b/ Hz1/2 - An obvious way 



-,1/2 



1/2 



to further decrease S^' and S,/ is to lower the SQUID 
inductance L, which can be done easily by decreasing 
the lateral distance between the JJs and by reducing the 
thickness of the Si0 2 layer separating the top and bot- 
tom Nb layers. In addition, the width of the Nb lines 
can be reduced further to increase (j)^ and to extend the 
range of magnetic fields where the SQUID can be oper- 
ated without vortices entering the wiring. All in all, we 
consider a spin sensitivity down to a few /.jb/Hz 1 ' 2 , for 
a field range exceeding 100 mT, to be achievable for this 
type of device. 



FIG. 4: (Color online) (a) I C (B) data of SQUID 2 for field 
sweep sequence 1-3 (black solid line) and 3-4 (red line plus 
symbols). Inset: Sj W (B) at optimum working point; dashed 
line is a linear fit. (b) Spectral density of rms flux noise 
Sl /2 (f) for B = OmT, 50 mT and 500 mT. 
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